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Abstract. The ratio between far-ultraviolet (FUV) and infrared (IR) luminosity densities from z — to z — 1 is 
discussed by using th e luminosity function s (LFs) of both wavelengths. The FUV LF [z = 0-1) based on GALEX 
has been reported bv lArnouts et alJ (|2005l) . whilst for the IR LF, we used t he IR AS PSCz 60-/im LF for the local 
universe and the Spitzer 15-pm LF at higher-z as used bv lLe Floc'h et al] (120051) . Both luminosity densities show 
a significant evolutionary trend, but the IR evolves much faster than the FUV. Consequently, the ratio pdust/pFuv 
increases toward higher-z, from ~ 4 (local) to ~ 15 (z ~ 1). It is also shown that more than 70 % of the star 
formation activity in the universe is obscured by dust at 0.5 z 1.2. 
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1. Introduction 

Dust attenuation is one of the most fundamental obsta- 
cles when we study the star formation activity of galaxies. 
Since ultraviolet (UV) radiation is emitted by young stars, 
it is in principle directly related to the recent star forma- 
tion rate (SFR). However, the use of UV light to trace the 
SFR is strongly hampered by the presence of dust, which 
absorbs and scatters the UV photons and finally re-emits 
the energy in the IR (mainly far-IR: FIR). Therefore, the 
UV and IR emissions play complementary roles in esti- 
mating the recent SFR of galaxies. 

In particular, the effect of dust has given rise to a 
long l asting debate o n the cosmic star formation history 
(e.g., iHonkinsl Eiol. The most direct way to address 
this issue is to compare the observed cosmic luminosity 
densities in UV and IR. In such studies, the luminos- 
ity function (LF) of galaxies is th e starting point (e.g., 
iTakeuchi. Yoshikawa. fc Ishiil EoOO). In this work, we in- 
vestigate the evolution of FUV and IR LFs. 

A highly reliable LF in the FUV has recently 
been published ba sed on new UV d ata obtained 
by GALEX 1 (e.g. lArnouts et alJ 120051) . In the IR, 
ITakeuchi. Yoshikawa. fc Ishiil ll2003bl) constructed a 
60 ^m LF from IRAS PSCz llSaunders et all 1200(1) . 
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Hereafter, we indicate IR emission integrated over 
8-1000 urn. by a subscript 'dust'. For higher-z, first 
results from S pi tzer 2 h ave been recently re porte d 
(iLe Floc'h et all l200Rt IPerez-Gonzalez et all 1200.1) . 
Making use of these LFs, we calculate the luminosity 
density at both wavelengths to examine the SFR history 
in the universe at < z < 1. 

Throughout this manuscript, we adopt a flat lambda- 
dominated cosmology with (h,Q ,\ ) = (0.7,0.3,0.7), 
where h = i?o/10Q[kms _ Mpc~ ], Slo is the density pa- 
rameter, and Aq is the normalized cosmological constant. 



2. Luminosity functions 

We define the luminosity function as a number density 
of galaxies whose luminosity lies between a logarithmic 
interval [logL,logL + dlogL]: 3 4>(L) = dn/dlogL. Here 
we define the luminosity at a certain wavelength band by 
L = vL,,. 



2.1. The FUV luminosity function 

IWvder et alJ l)2005|) estimated the local LF of galaxies at 
FUV (1530 A) from GALEX data in overlapped regions 
with 2dFGRS l|Colless et alJl200l|) . The local FUV LF is 
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Table 1. Schechter parameters for FUV luminosity func- 
tion (converted from lArnouts et afll2005j) . 



Redshift 


Q 


L„(FUV) 
h- 2 [L ] 


0. 

h 3 [Mpc" 







1.21 


1.81 x 10 y 


1.35 x 10 


-•i 


0.2-0.4 


1.19 


2.43 x 10 9 


1.95 x 10 


-2 


0.4-0.6 


1.55 


6.75 x 10 9 


5.38 x 10 


-3 


0.6-0.8 


1.60 


9.32 x 10 9 


5.25 x 10 


-3 


0.8-1.2 


1.63 


1.19 x 10 10 


3.63 x 10 


-3 



well described by a Schechter function l|Schechterlll976|) 



4>{L) = (In 10) 



exp 



(1) 



(1.21,1.81 x 10 9 ^ 2 L0,1.35 x 



At z = 0, (a, £»,</>») 
lQ - 2 h 3 Mpc~ 3 ). 

lArnouts et alJ l|2005|) presented the evolution of the 
GALEX FU V LF using the VIRM OS VLT Deep Survey 
(VVDS: see, iLe Fevre et ail l200.il . They found that the 
FUV LFs at hig her z ar e also well fitted by the Schechter 
function. lArnouts et al.l l)2005|) directly measured the pa- 
rameters (a, .L*, </>*) at each redshift bin. They reported 
that the a and L» monotonically increase with z, while 
the (ft* decreases with z. We adopt the converted values of 
the parameters in Equation (Table EJ. 

2.2. The IR luminosity function 

Contrary to the FUV LF, the local IR LF is not well-fitted 
by a Schechter fu nction, although it can be expressed as a 
function given bv lSaunders et all l)l990fl which is defined 

as 



0(L) = (f>* I — 



exp 



1 

2a 2 



1 2 



(2) 



We use the parameters for t he loca l 60 /im LF, given by 
iTakeuchi. Yoshikawa. fc Ishiil l|2003b|) . We converted the 
60 /im luminosity to that of the total dust emission by 
adopting the average 60/im-to-d ust flux ratio of 2.5 esti - 
mated from the PSCz sample of ITakeuchi et all l|2005a|) . 
This luminosity-independent conversion can be justified 
because of the tight linear correlation of the two quanti- 
ties (correlation coefficient r — 0.991). 

Since most of the galaxies in IRAS PSCz are local 
(z < 0.1), we need a deep survey result to evaluate the 
evolution of the IR LF. Recently, very important results 
from Spitzcr MIPS 24- /im observations ha ve been reported 
on the mid-IR (MIR: 12 or 15 u rn) LF l|Le Floc'h et alJ 
2005: IPere z- Gonzalez et al .120051) . These authors reported 
a very str ong evolution tr e nd for the MIR LF. Hereafter 
we adopt ILe Floc'h et alJ l)2005|) because their LFs are 
given in a form identi cal to the one adopted by us. 
However, we remark that lPerez-Gonzalez et al . ( 2003), al- 
though adopting a diff erent form to our study and that of 



ILe Floc'h et al. l)2005|) . nonetheless also reached consis- 



iLe Floc'h et alJ |2005) nrs ^ estimated nonparametric 
15-/xm LFs for each redshift bin at < z < 1 from MIPS 
24- /im data by the 1/Vmax method. Then, using some 
spectral energy distribution (SED) templates, they con- 
verted their 15-/iin luminosity to the total dust luminos- 
ity. Based on these nonparametric LFs at < z < 1, they 
estimated the evolution rate of the IR LF, by adopting 
the form 



4>{L,z) = g(z)<j) 



(3) 



where (j>o(L) is the local functional form of the LF. They 
assumed a power-law form for the evolution functions f(z) 
and g(z) as 



f(z) = (l + zf, g(z) = (l + zY 



(4) 



Through a % 2 minimization between the nonparametric 
LFs and Equation they obtained the parameter values 



h0.7 



and P = 0.7^n R- This means that £* 



oc (1 + z) - 7 in Equation (J2J, whilst a 



cis Q — 3.2 q 2 

(1 + z) 3 - 2 and 
remains constant. 

Here we mention the uncertainty in the conversion of 
ILe Floc'h et al.l 1120051) converted th e local 
Takeuchi. Yoshikawa. fc Ishiil l|2003bj) to the 
dust LF by using their model SED templates. Since their 
conversion procedure is different from ours, we examined 
their consistency. The difference between their conversion 
and ours does not exceed 5-%, thus we judge the difference 
to be negligible for the subsequent analysis at 60 /im for 
z = 0. 

At higher-z, they convert L15 to Ldust by SED tem- 
plates. We also consider the potential systematic un- 
certainty introduced by this procedure. In particular, 
the evolution of the population causes an increase of 
the fraction of IR luminous galaxies (LIRGs, ULIRGs), 
which may have different SEDs to less active galaxies. 
This can lead a systematic change of the correspond- 
ing template SEP wit h z. To evaluate this uncertainty, 
ILe Floc'h et all §005) made a comparison between sev- 
eral IR SED templa t es llDale et al.ll200lHCharv fc Elbazj 
200l | iDale fc Helm J El02t IChaniall l200.j lLagache et alJ 
20031 12004|) . Fromtheir Fig. 8, they estimated the typ- 
ical uncertainty to be ~ 0.2 dex. We further extended 



their test using some other SED librar ies ( Efstathiov^et^l 
2000tlEfstathiou fc Rowan-Robinsonl2003t ITakeuchi et al 
latt 



tent conclusions on the amount of the evolution. 



2001albh . and confirm their claim. The good linearity be 
tween MIR (IRAS 12 and 25 /im) and dust luminosities 
at a very wide range of luminosity (10 6 L Q -10 13 L Q ) 
l|Takeuchi et alJ l2005al). even for galaxies wit h extreme 
SEDs (see, e.g.. iTa.keuchi et a.lJl2003al I2005W) . also en- 
sures the robustness of the estimation. In summary, the 
uncertainty of the dust luminosity estimation is a factor 
of three. We should keep this uncertainty in mind for the 
following. 
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Fig. 1. The evolution of the luminosity function (LF) of 
galaxies in the far-ultraviolet (FUV: 1530 A) obtained by 
GALEX and infrared (dust: 8-1000 fim) obtained from 
IRAS PSCz (z = 0) and Spitzcr at higher-z. Thick lines 
show the LFs in the FUV, and thin lines depict those of 
dust. Vertical tick marks on the LF indicate the lowest 
luminosity limits above which the observed data exist. 
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Fig. 2. The evolution of the contribution of galaxies to 
the luminosity density in the FUV and IR. The meaning 
of the different lines are the same as that of Figure H 



3. Result and discussion 

3.1. Evolution of the luminosity contribution in the 
FUV and IR 

The evolution of the FUV and dust LFs are shown in 
Figure ^ First, it is worth mentioning the (well-known) 
difference of the local LF shape of FUV and dust: for the 
dust LF, bright galaxies (L ^ 10 10 L©) are much more 
numerous than those in the FUV. This leads to the dif- 
ference in the main population contributing to the total 
emitted energy. The product Lip(L) represents the energy 
contribution of galaxies with luminosity L (see Figure |5J) . 
In the FUV, the main contributor is L» galaxies, with 
fainter galaxies emitting a non-negligible fraction of en- 
ergy at z > 0.5. In contrast, the effect of the evolution ap- 
pears in the bright end for the dust LF. The contribution 
from the most luminous galaxies increases with redshift. 

3.2. Evolution of the FUV and dust luminosity 
densities and the mean dust attenuation 

For FUV, we integrate L<j){L) over Lfuv = 
10 6 L Q -10 15 Lq at each z to obtain the evolution 
of the FUV luminosity density pfuv, while for the 
dust, we first construct a LF at a given z according 
to Equations ( I T?f> n and ijijl with the estimated value of 
iLe Floc'h et alJ l(2005h . and integrate L<j)(L) over the 
same range as that of FUV galaxies. The densities depend 
very little on the integration luminosity range: even if we 
change the lower bound to 1 Lq, the integrated value 
only increases by 0.4 % for the dust luminosity and by 
3 % in the FUV luminosity. 

The luminosity densities are summarized in Table [21 
Both luminosity densities show a significant evolution- 
ary trend, but the dust luminosity density evolves much 
faster than that of the FUV. In Table El we only 
list the systematic uncertainty potentially introduced 
by the c hoice of SEP templates. For statistical er- 
rors, see Schiminovic h et al ] l)2005|) for the FUV and 
ILe Floc'h et alJ l)2005h for 15 /zm. Consequently, the ra- 
tio Pdust/pFUV increases toward higher-z, from ~ 4 (local) 
to ~ 15 (z ~ 1), i.e., the dust luminosity dominates the 
universe at z ~ 1. 

The FUV-to-dust luminosity density ratio can be in- 
terpr eted as the global m ean dust attenuation in the uni- 
verse. iBuat, et, al.l (j2005li provided a formula which relates 
the dust to FUV flux ratio to the dust attenuation in the 
FUV, A(FUV) [mag], as 

A(FUV) = -0.0333?/ 3 + 0.3522y 2 + 1.1960?/ + 0.4967 , (5) 



where y = log (Fd us t/FFuv) and F — vS v (S u - flux den- 
sity) . The mean attenuation obtained by Equation (J5J is 
also tabulated in Tabled 
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Table 2. Evolution of the FUV and dust luminosity densities, the mean dust attenuation, and the converted cosmic 
SFR densities. 



Redshift p FU v Pdust Pdust/pFW A(FUV) psfr(FUV) p S FR.(dust) Hidden SFR 

h [L Mpc~ 3 ] h [L©Mpc~ 3 ] [mag] h [Moyr^Mpc -3 ] h [Afoyr^Mpc -3 ] fraction 
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3.3. Fraction of obscured star formation 

We interpret the data in terms of SFR. Assuming a con- 
stant SFR o yer 10 s yr, and Salpeter initial mass func- 
tion (IMF) jSaroeterl Il955l mass range: 0.1-100 M Q ), 
Starburst99 ( Leitherer et aljfl999j) gives the relation be- 
tween the SFR and L(FUV) = vL v at FUV (1530 A), 

log i(FUV) = 9.51 + log SFR . (6) 

For the IR, to transform the dust emission to the SFR, 
we assume that all the stellar light is absorbed by dust. 
Then, we obtain the following formula under the same 
assumption for both the SFR history and the IMF as those 
of the FUV, 

log i(dust) = 9.75 + log SFR . (7) 

However, a significant fraction of the dust emission is due 
to the heating of grains by old stars which i s not d irectly 
related to the recent SFR. iHirashita et alJ (|2003j) found 
that about 40 % of the dust heating in the nearby galax- 
ies comes from stars older than 10 s yr. Adopting this 
correction, we obtained the evolution of the star forma- 
tion rate densities from FUV and dust (/osfr(FUV) and 
PSFR(dust)) which are presented in Figure|2|and Tabled 
The evolution of pdust and pfuv (therefore psFR.(dust) 
and Psfr(FUV)) is described by the form of (l + z) R . The 
power-law index R of pdust is estimate d to be 3.9 ± 0.4 by 
iLe Floc'h etaH l|2005|) . and for 

Pfuv > ISchiminovich et all 
l)2005|) give an index of 2. 5 ±0.7 (dotted and solid curves in 
Figure (5J. Th e evolution of pdnst is slightly stronger than 
suggested by iLaeache et all 1120031) . but consistent with 

" l200lh and lTakeuchi et alJ 
2005) adopted a linear re- 



those 



given bvlCharv fc Elbazf i 
<|2001afl . lPerez-r T ouzalez et TIT I 



lation between th e SFR and the dust luminosity given by 
iKennicuttl l[l99St) and obtained the evolution of the SFR 
density as /osFR(dust) oc (1 + z) 4 0±0 ' 2 , which is in perfect 
agreement with the above result for the dust luminosity. 
Therefore, this is a robust conclusion which does not de- 
pend on the adopted LF functional shape for the fitting 
or the details of the conversion from dust emission to the 
SFR. 

The fraction of obscured SFR density is also presented 
in Table [5] About half of the star formation activity is 
obscured in the local universe, while at 0.5 ^ z ^ 1.2, 
about 70 % of the total SFR is hidden by dust. In partic- 
ular, the obscured SFR fraction reaches more than 80 % 
at z ~ 1. The result is consistent with a previous sug- 
gestion from the optical-to-FIR luminosity density ratio 
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Fig. 3. The evolution of the star formation densities de- 
rived from FUV and dust. Typical statistical error is about 
a factor of two, and systematic uncertainty related to IR 
SED template is about ~ 0.2 dex. 



(e.g., iTakeuchi et"al"ll2001cj) . This result should substan- 
tially change the way we see the SFR history at high-z. 
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